Static potential and local color fields in unquenched lattice QCD3 
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String breaking by dynamical quarks in three-dimensional lattice QCD is analyzed through mea- 
surements of the potential and the local color-electric field strength generated by a static quark- 
antiquark pair. Simulations were done for unquenched SU(2) color with two flavors of staggered 
light quarks. An improved gluon action was used, which allows simulations to be done on coarse 
lattices, providing an extremely efficient means to access the large quark separations and long prop- 
agation times at which string breaking occurs. The static sources were generated using Wilson 
loop operators, hence no light valence quarks are present in the resulting trial states. Results give 
unambiguous evidence of string breaking. First the static potential is shown to saturate at twice the 
heavy-light meson mass at large separations. Then it is demonstrated that the local color-electric 
field strength in the region between the heavy quarks tends towards vacuum values at large sepa- 
rations, the first time that this most graphic effect of quark vacuum polarization on the confining 
flux-tube has been realized in lattice QCD. Implications of these results for unquenched simulations 
of four-dimensional QCD are drawn. 



I. INTRODUCTION 

The string model is widely used in phenomenologies 
of quark confinement and hadronization. A basic sys- 
tem of interest in the string model consists of a heavy 
quark (Q) and antiquark (Q) pair, which are assumed 
to be bound by a narrow tube of color flux. This model 
naturally gives rise to a linearly confining potential, if 
the tube has a cross-section that is approximately inde- 
pendent of the quark separation In fact simulations 
of quenched lattice QCD long ago demonstrated the for- 
mation of color-fiux tubes having properties that are in 
accord with simple models. 

In full QCD the confining QQ string should not per- 
sist to arbitrarily large separations R, due to the effects 
of dynamical sea quarks. One expects that it becomes en- 
ergetically favorable for a light quark (g) and antiquark 
(g) to materialize from the Dirac sea at large R, and the 
light quarks should bind to the heavy quarks resulting 
in the formation of a pair of noninteracting color-neutral 
mesons. This phenomenon of hadronization, or string 
breaking, implies for instance that the potential V{R) 
for a static QQ pair should saturate at large separations. 



Hm V{R) = 2Mq^, 

R^oo 



(1) 



where Mqq is the mass of the ground-state heavy-light 
meson. 

Observation of Eq. |^ in an unquenched lattice QCD 
simulation provides the simplest criterion to demonstrate 
string breaking. A much more detailed probe of this basic 
phenomenon of hadronic physics would come from mea- 
surements of the local color fields f^iy{x; R) as functions 
of position x in the region between the static valence 
quarks, for various values of R. The effects of dynami- 
cal quarks should be to suppress the fields in the region 
between the static sources, with the fields saturating at 
vacuum values in the limit of large R 



lim f^y{x;R) = 0, 



(2) 



for X outside a finite region surrounding each static va- 
lence quark, corresponding to the size of the Qq meson. 

Despite considerable effort, however, it has proven to 
be very difficult to estabfish Eq. J^l, and especially Eq. 

in unquenched lattice QCD simulations 2,'^,'4|. The 
traditional computational strategy is to use Wilson loop 
operators in order to create trial states containing only 
static sources, with no valence light quarks. In that case 
Eqs. (P) and (O provide unambiguous evidence of the 
effects of light quarks that "materialize" from the sea 
(one must take care to note however that Euclidean lat- 
tice simulations only provide an analogy to the real-time 
process of hadronization) . The difficulty in making these 
measurements on the lattice is due in part to the expo- 
nential suppression of the correlators, reflecting the large 
energy of the state at large separations R, and the large 
propagation time T which is necessary in order to isolate 
the ground state of the system 0, 0, H ■ 

A convincing demonstration of the saturation of the 
unquenched V{R) has only recently been achieved, 
thanks in part to the realization that improved actions 
can be used to do accurate simulations on lattices with 
coarse grids, where the computational effort can go into 
generating large ensembles of gauge-fields at the large 
length scales relevant to string breaking, rather than to 
generating short-distance modes that are not of interest 
here |^ U Studies of the local color-fields are still more 
problematic than the static potential; the simplest defini- 
tion of ff_ii,{x; R) comes from the correlator between the 
Wilson loop and the square of the local field strength, 
which requires a vacuum subtraction, leading to a very 
poor signal-to-noise ratio. A demonstration that the 
fields in the region between the heavy quarks become 
vacuum saturated has not been realized to date. 

In this work we turn to Euclidean lattice QCD in three 
dimensions (QCD3), exploiting the computational ad- 
vantage of this toy model (including a power law sup- 
pression of lattice discretization errors in this super- 
renormalizable theory) in order to convincingly demon- 
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strate that dynamical quarks lead to saturation of the 
static potential and the local color fields. This has been 
done using very little computational power. 

Lattice studies of confinement physics using QCD3 
have a long history; this theory exhibits most of the 
hallmarks of four-dimcnsional QCD, including confine- 
ment, flux-tube formation, chiral symmetry breaking, a 
rich glueball spectrum, and a finite-temperature decon- 
fining phase transition 

Some of the work presented here is a follow-up to an 
earlier study on string breaking in unquenched QCD3 by 
one of us 0, Q . What is new includes a much more ex- 
tensive study of the static potential on larger lattices and 
with higher statistics, enabling measurements of the po- 
tential out to substantially larger R. More important, 
in this work we provide the first demonstration in un- 
quenched lattice QCD, as defined in any number of di- 
mensions, of vacuum saturation of the local color fields. 

Although these results are for a three-dimensional the- 
ory, they have important implications for string breaking 
studies in QCD4. In particular we find that string break- 
ing is observed using trial states without light valence 
quarks, when they are propagated for Euclidean times 
T w 1 fm, longer than what has been attained in earlier 
unsuccessful studies, but not excessively so, this being the 
natural scale associated with hadronic binding. [As dis- 
cussed in Sect.^we find a consistent relation between the 
physical length scales in QCD3 and in four-dimensional 
hadronic physics, using two very different quantities to 
identify the lattice spacing in QCD3: the quenched string 
tension, and the lightest unquenched vector meson mass, 
analogous to the p.] 

Saturation in both the potential (as previously ob- 
served in 0), and in the local color fields, is found to 
occur at separations R ^ 1.5 fm, which is similar to es- 
timates of the string breaking scale in QCD4 The 
correlation functions are noisy, yet a sufficiently large en- 
semble of configurations was generated with relatively lit- 



tle computational effort, even by the standards of the in- 
trinsically cheaper three-dimensional theory, because we 
used improved actions to accurately simulate on coarse 
lattices, with spatial spacings as ~ 0.2 fm. The same 
reasoning should apply to the situation in unquenched 
QCD4 6]. 

The rest of this paper is organized as follows. In Sect, 
^we describe the lattice actions and simulation param- 
eters, and we also discuss our scale setting procedure for 
making contact between QCD3 and QCD4. Results for 
the static potential are presented in Sect. UTTI and results 
for the local field strength in Sect. II Vl A summary and 
some further discussion of the implications of this work 
for simulations of QCD4 are found in Sect. IVl 



II. SIMULATIONS 

Unquenched simulations of SU(2) color in three dimen- 
sions were done using a tree- level 0(a^ )-accurate im- 
proved gluon action, allowing for different "temporal" 
and "spatial" lattice spacings at and as, respectively 



J2 

X.fJ,>V 



3 



12 



(3) 



where P^^ is one-half the trace of the 1x1 plaque- 
tte and R^y is one-half the trace of the 1x2 rectan- 
gle in the fi x v plane. The bare lattice anisotropy is 
input through ^.^ = £,q = (at/as)bare [for «,j = 1,2], 
and ^fc3 = ^3fc = 1/^0- The three-dimensional theory is 
super-renormalizable, with a bare coupling gf, having di- 
mensions of mass, and which enters Eq. through the 
dimensionless combination (3 = A/{g'^as)- 

We use the unimproved Kogut-Susskind staggered- 
quark action, which in three dimensions describes two 
flavors of four-component spinors illij 



J 



S'ks = 2amo ^ x{x)x{x) -f ^ C,t,-ntJ,{x)x{x) [Uf,{x)x{x + (i) 



Ulix- 



fi)x{x - A)] , 



(4) 



where ry^(x) = (— 1)^i+---+^m 
phase, and where (^1,2 = 1, (3 



^ is the usual staggered 
1/Co. 



Simulations were done on a 22^ X 28 lattice at /3 = 3, 
with a bare anisotropy ^0 = 1/2, and a bare quark mass 
TOo/(7o =0.1. We find that this corresponds to a rather 
large pion mass, with rrij^/mp ~ 0.61. Hence the light 
quark used here is actually rather heavy, comparable to 
the strange quark; nevertheless this quark mass is ap- 
parently light enough to resolve string breaking at acces- 
sible separations R. These parameters are very similar 
to those that were used in Ref. although the lattice 



volume is substantially larger here. For comparison the 
simulation parameters in the present study are compared 
with those in Ref. ||] in Table |T| 

The configurations were generated using the hybrid 
molecular dynamics algorithm (the ^-algorithm) 
with time step At — 0.02 and 50 molecular dynamics 
steps taken for each trajectory. Direct calculation of au- 
tocorrelation times for the largest Wilson loops showed 
T < 0.2 and, based on binning studies of autocorrelations, 
we elected to skip 10 trajectories between measurements. 
Our final data set consists of 30,000 measurements, cor- 
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Simulation Parameters 


Ref. f5] 


This work 


Lattice Volume 


\& X 10 


22^ X 28 


P 


3.0 


3.0 


(fltj t[s)barc 


isotropic 


1/2 




0.075 


0.10 


incas 


6,000 


30,000 



TABLE L Comparison of simulation parameters used in an 
earlier study by one of us Q, and in the present work. The 
quenched and unquenched simulations in this work were done 
at the same couphngs. A'^moas is the number of measurements 
in each unquenched study. 



responding to a total of 300,000 trajectories. 

The bare parameters go, "^o, f^nd are finite in the 
continuum limit of this super-renormalizablc theory. At 
finite lattice spacings these parameters must absorb fi- 
nite renormalizations in order to keep physical quantities 
fixed. However in Refs. |^ ^3 it was found that little 
renormalization is required when the improved gluon ac- 
tion is employed, in the range of couplings used here. 
For example, the unquenched static potential showed no 
discernable change from /3 = 3 to /3 = 2, corresponding 
to a 50% increase in lattice spacing 0; meson masses 
were also found to have small discretization errors |l2j |. 
Of relevance to our comparison of the quenched and un- 
quenched results is the fact that the lattice spacings in 
the two theories are found to be very similar at a given 
bare coupling fi (a quantitative comparison is made be- 
low). Moreover one finds similar meson masses in the 
two theories, at a given bare quark mass; for example we 
find m-j^lmp ~ 0.75 in the quenched theory, compared to 
the value w 0.61 unquenched. 

We likewise find very little renormalization of the bare 
lattice anisotropy. The renormalized anisotropy was de- 
termined by comparing measurements of the static po- 
tential where the direction of Euclidean time was oriented 
along the at lattice axis and then along an axis. The 
results of the quenched simulations give a renormalized 
value at I as = 0.478(7), while in the unquenched theory 
we find at I as — 0.455(9) (much larger renormalizations 
are found if the unimproved Wilson gluon action is used 



i). 



One can make rough contact between the length scales 
relevant to string breaking in QCD3, and the corre- 
sponding scales in four-dimensional QCD, by identify- 
ing a confining scale in the two theories For in- 
stance, one can identify the quenched string tension in 
QCD4, ^/a ~ 0.44 GeV, with the result in quenched 
three-dimensional SU(2) color [1], ^ j gl = 0.3353(18). 
In effect this allows one to express the dimensionful cou- 
pling l/f^Q, or equivalently the lattice spacing at a 
given /3, in the three-dimensional theory in "fermi." For 
instance at /3 = 3 one finds as ~ 0.20 fm. The con- 
sistency of this procedure can be established by iden- 
tifying other physical quantities in QCD3 and QCD4. 
We measured the mass my of the lightest vector meson, 



with the result atmy = 0.425(8) in the unquenched case, 
and identifying this with the physical p meson mass one 
finds tts ~ 0.24 fm. In the quenched theory at the same 
couplings the vector meson mass is atmy — 0.380(2), 
which would identify the spacing as as ~ 0.20 fm. [This 
also supports the contention that the quenched and un- 
quenched QCD3 theories can be directly compared at 
the same bare couplings, though a small mismatch in the 
couplings will not affect the trends that we will use to 
establish Eqs. (QJ and ^] 

It is important to note at this point that one can, and 
we will, establish the string breaking criteria completely 
within the QCD3 theory, independent of this identifica- 
tion of scales with QCD4. This connection is only used to 
the extent that it gives one some physical intuition for the 
length scales that are being probed in QCD3, and since 
one may also thereby make some comparison between the 
parameters of the simulations that are done here, which 
successfully resolve string breaking, and those that have 
been used in various studies in QCD4. It is to these 
ends that we sometimes label our simulation parameters 
according their values in "fermi." We return to the com- 
parison of our simulation results with QCD4 in Sect. Ivl 



III. RESULTS: STATIC POTENTIAL 

We extracted the static potential from expectation val- 
ues {W{R,T)) oi RxT Wilson loops. Standard APE- 
type smearing was used for the spatial links in the Wilson 
loops, in order to enhance the overlap with the ground 
state. Measurements were done both for on-axis and 
many off-axis separations (i.e. using link paths joining 
lattice sites displaced along both spatial directions). 

In order to demonstrate the onset of string breaking it 
is essential to carefully study the time-dependent effec- 
tive potential V{R, T) as a function of the propagation 
time T, where 



V{R,T) = - In 



{W{R.T)) 

(i¥(i?,r-i)) 



(5) 



with the true ground state energy V{R) given as usual by 
V{R) = limT^oo V{R, T). Results for the unquenched ef- 
fective potential are presented in Fig.^for several values 
of T, and are compared with the quenched potential (the 
latter is evaluated at T = 12at, which is generally large 
enough to isolate the ground state in that case, except at 
the largest i?'s in the figure). The potential is also com- 
pared with the unquenched simulation determination of 
twice the heavy-light meson mass (using an appropriate 
staggered interpolating field 14]). 

These results clearly demonstrate string breaking as 
defined by Eq. QJ. An asymptote in the unquenched 
potential at large separations is progressively revealed 
as T increases; at the longest propagation time shown 
(T = 12at w 1.2 fm), the potential has clearly saturated 
at the expected limit of 2MQq over a very wide range of 
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FIG. 1: Time-dependent unquenched effective potential 
V{R, T) versus quark separation R, at several propagation 
times T. The quenched potential is also shown. The dashed 
lines give the unquenched simulation measurement of 2MQq. 



radii. From these results one can make a rough determi- 
nation of the separation i?sb for the onset of string break- 
ing, corresponding to the point at which the potential 
reaches the correct asymptote. We find Rsh ~ (6-8)as. 
Using the scale setting procedure described in Sect. ^ 
to convert the lattice spacing to "fermi," one obtains 
^sb ~ (1.2-1.6) fm, which is very similar to estimates 
of the string breaking distance in QCD4 ^loi . 
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FIG. 2: Time-dependent unquenched 
sion a^a{T), dehned according to Eq. 
show the quenched string tension. 




"effective" string ten- 
@. The dashed lines 



To further analyze the approach to an asymptote in 
the effective potential as T is increased, we define an 
"effective" string tension cr(T), computed from V{R,T) 
over a small interval in separation R that is beyond the 
string breaking distance Rsh] we take 



Results for a{T) so defined were obtained from fits to 
the simulation measurements, and are shown in Fig. |21 
The unquenched string tension clearly tends towards zero 
with increasing T, and this trend is manifest even at 
the smallest propagation times, where the measurements 
are most accurate. An extrapolation by-eye of the trend 
in a{T) from the smallest times would suggest that the 
string tension vanishes at a propagation time that is only 
slightly longer than what is revealed by direct inspection 
of an asymptote in V{R, T) at large T, which clearly sets 
in by T = 12at, as seen in Fig. ^ 



IV. RESULTS: LOCAL COLOR FIELDS 

The simplest observable for measuring the local color 
fields comes from the correlator of a Wilson loop W{R, T) 
with the plaquette Pfiu 



{W(R,T)P^,{x)) 
{W{R,T)) 



{P. 



fJ.V / 



(7) 



where we suppress the dependence of f^i, on R and T. 
We measure x relative to the center of the Wilson loop. 
In the continuum limit (and at sufficiently large T) the 
above correlator corresponds to the expectation value of 
the square of the Euclidean field strength in the presence 
of the static QQ pair, after vacuum subtraction 

\im UA^)^(TT{F^,ix)f) -hi-{F^,f) , (8) 

\ / QQ \ / vac 

where F^,^ denotes the usual continuum field strength, 
and (. . ■)qq denotes the expectation value in the pres- 
ence of the static sources, and (. . .)vac is the vacuum ex- 
pectation value. 

The most important component of the field correlator 
corresponds to the color-electric field with vector com- 
ponent along the direction parallel to the axis joining 
the quark and antiquark (which, for instance, makes the 
largest contribution to the potential energy) . We denote 
this correlator by E^^(x), given by 



E\\ix)^~hti£), 



(9) 



V{R,T) ^a{T)R + b, i?/a, e [8, 10] 



(6) 



where the minus sign accounts for the Euclidean contri- 
bution to the energy density, and where here the "1" di- 
rection denotes the QQ axis and t denotes the Euclidean 
time axis corresponding to the small lattice spacing at- 

We present one-dimensional profiles of the field 
strength, along two different lines through the flux-tube, 
in Figs. island^ We write x = (a;||,a;_L), where a;|| is the 
position along the central axis joining the quarks, and x± 
is the position perpendicular to that line, both measured 
from the center of the axis. We plot the profile along the 
central axis joining the quarks in Fig.O and along a line 
perpendicular to the central axis in Fig. 01 The quenched 
and unquenched correlators are shown in each case and, 
in order to emphasize the shape of the profiles in the two 
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FIG. 3: Color-electric field _B"(a;||,a;x = 0) as a function of FIG. 4: Color-electric field -E"(a;|| — 0,x±_) as a function of 



position x\\ along the axis joining the quark and antiquark, 
for two values of the separation, at a fixed T. -Bmax is the 



position x±_ along a line perpendicular to the axis joining the 
quark and antiquark, for two values of the separation, at a 



maximum field strength, in the given theory, along this line fixed T. iS^ax is the maximum field strength, in the given 



(i.e. the field at the location of either quark). 



theory, along this line (which in this case corresponds to the 
field at the midpoint of the QQ axis) . 



theories, each correlator is normahzed to its maximum 
value along the particular line. 

The existence of a flux-tube is demonstrated in part 
by a "plateau" in the field strength in a region between 
the sources (which is expected to increase in width as R 
increases), where the field is approximately independent 
of xy, as seen in Fig. O The relative strength of the field 
in this region is suppressed in the unquenched theory. 
[We note that i?" in the plaquettc nearest to a static 
source, at |a;||| = int(i?/2), x±_ — 0, is generally about 
half as large in the unquenched theory as in the quenched 
case, at large R. To the extent that the couplings in the 
two theories are well matched, this already illustrates the 
suppression of the fields due to the effects of unquenching. 



since the correlator is an average of the squared-field over 
a unit cell.] 

Notice as well the approximate symmetry of the fiux- 
tube profiles about the center of the QQ axis. The 
transverse widths of the flux tubes can be inferred from 
Fig. ^ and are indistinguishable in the unquenched and 
quenched cases (the width is also roughly independent of 
separation, which is consistent with simple models 1]). 

The results in Figs. 13 and 21 are shown for fixed T = 
7 at ~ 0.7 fm. As indicated by the analysis of the poten- 
tial in the previous section, one can easily be misled as to 
the extent of string breaking unless one attains propaga- 
tion times of at least 1 fm. This is particularly difficult 
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in the case of the field correlators. Alternatively one can 
measure at several "intermediate" values of T, and study 
the systematics of increasing propagation times. To il- 
lustrate this we compare the fields at the midpoint of the 
QQ axis in the two theories 

p{R,T)=\e\\{x = 0)] /\e\\{x = Q)\ 

L J unquenchcd / L J quenched 

(10) 

where we explicitly indicate the dependence of the ratio 
on R and T . The results are shown in Figs.[Sland|Hl which 
provide analogues for the fields to the effective potential 
plots in Figs, n] and 121 




FIG. 5: Ratio of the field strengths, unquenched to quenched, 
at the midpoint of the QQ axis, as a function of separation 
R, for various propagation times T. 




FIG. 6: Ratio of the field strengths, unquenched to quenched, 
at the midpoint of the QQ axis, as a function of propagation 
time r, for fixed R = 8as. 

The unquenched fields are suppressed at all R and T, 
however at modest propagation times the unquenched 



field does not drop with increasing R; taken at face value, 
this might suggest a failure to observe string breaking, 
perhaps due to some pathology in the trial states gen- 
erated by Wilson loops, as has been suggested by other 
authors (see e.g. Refs. 0,01). In actuality this is sim- 
ply a reflection of the fact that the ground state is not 
resolved unless one achieves propagation times « 1 fm, 
which is characteristic of hadronic binding. 

The systematics are particularly well illustrated in Fig. 
El which shows the time-dependence of the field strength 
for a fixed separation R — Sa^, which was chosen because 
it is just beyond the string breaking distance Rsh- The 
trend is once again manifest even at small propagation 
times, where the signal is very clean, and an extrapola- 
tion by-eye indicates that the field becomes vacuum satu- 
rated for T > 12at, which is consistent with our analysis 
of the effective potential in the previous section. 



V. SUMMARY AND OUTLOOK 

We were able to observe string breaking by dynamical 
quarks in three-dimensional QCD by working on a coarse 
lattice, using the staggered quark action along with an 
improved gluon action. String breaking was observed at 
large quark-antiquark separations, both as saturation of 
the static potential at twice the heavy-light meson mass, 
according to Eq. (Q, and as saturation of the local field 
strength at vacuum values, according to Eq. Mea- 
surements were done using only trial states generated by 
Wilson loop operators, in which no valence light quarks 
are present. Hence saturation of the unquenched poten- 
tial and fields are highly nontrivial results, due entirely 
to the effects of dynamical sea quarks. The local field 
strength is particularly difficult to compute at the prop- 
agation times and separations at which string breaking 
occurs. Nonetheless computations were done over a range 
of propagation times that was sufficient to clearly estab- 
lish the trend for the field strength to approach vacuum 
values in the region between the static sources. 

Although this work was done in QCD3, the results pro- 
vide clear implications for large scale simulations of full 
QCD in four dimensions 5, 6]. Of particular importance 
is the need to achieve sufficiently long propagation times 
for the correlation functions; one can easily be misled as 
to the extent of string breaking unless one attains prop- 
agation times of at least 1 fm, the scale characteristic 
of hadronic binding. It is also sensible to work on rel- 
atively coarse lattice, using improved actions to achieve 
accurate results, since this alleviates the computational 
burden associated with generating short-distance modes 
that are not relevant to string breaking. 

In this connection we compare the parameters of two 
simulations of unquenched QCD4 (which used trial states 
generated only by Wilson loops) and of the present study 
of QCD3, in Table ITU These studies are distinguished in 
part by whether the string breaking condition Eq. was 
resolved. One must exercise caution when comparing the 
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QCD4 Ref. m] 


QCD4 Ref. ^ 


QCD3 (this work) 


Eq. Q? 


No 


Yes 


Yes 


as 


0.15 fm 


0.40 fm 


0.20 fm 




0.59 


0.25 


0.61 




0.4-0.8 fm 


1.2 fm 


1.2 fm 


-^max 


1.8 fm 


1.2 fm 


2.8 fm 


^ mcas 


492 


1,000 


30,000 


Volume 


16=^ X 32 


6^ 


22^ X 28 



TABLE II: Comparison of some recent unquenched lattice 
simulations in QCD4 and in the present QCD3 analysis, dis- 
tinguished in part by whether the particular study resolved 
string breaking in the static potential according to Eq. Q. 
The lattice spacings, meson masses, maximum quark separa- 
tions -Rmax and correlation function propagation times Tmax, 
number of measurements A'mcas , and lattice volumes are com- 
pared. 



two simulations in QCD4, which were done with different 
actions and simulation algorithms. One must likewise ex- 
ercise caution in the comparison between the QCD4 and 
QCD3 simulations (taking note of the discussion at the 
end of Sect . ITllconcerning the scale setting procedure that 
is used for this purpose). Nonetheless this comparison is 
instructive, as it highlights the advantages to be gained 
by working on coarser lattices. In unquenched QCD4 the 



computational cost scales with lattice spacing roughly as 
a^^, hence a small increase in the coarseness of the lattice 
can significantly reduce the computational cost, allowing 
one to generate much larger ensembles and/or to work 
at much smaller quark masses. Indeed the QCD4 study 
in Ref. 6] , which did resolve string breaking in the static 
potential, only required the equivalent of few PC-years of 
run time (comparable to the cost of the present study) . 

In this work we balanced the desire to use the coars- 
est lattice possible, in order to minimize the computa- 
tional cost, against the need to have a spacing that is 
fine enough to resolve the spatial distribution of the color 
fields. We also systematically studied the trends in cor- 
relation functions with increasing propagation times, in 
order to leverage the cleaner data at smaller T. It is 
reasonable to anticipate that string breaking in the lo- 
cal color fields, Eq. Q, may be accessible to simulations 
in QCD4, following the approach that proved successful 
here in QCD3. 
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